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ABSTRACT 

We use a cosmological numerical simulation to study the tidal features produced by a minor merger 
with an elliptical galaxy. We find that the simulated tidal features are quantitatively similar to 
the red tidal features, i.e., dry tidal features, recently found in deep images of elliptical galaxies at 
intermediate rcdshifts. The minor merger in our simulation does not trigger star formation due to 
active galactic nuclei heating. Therefore, both the tidal features and the host galaxy are red, i.e. a dry 
minor merger. The stellar mass of the infalling satellite galaxy is about 10 10 M , and the tidal debris 
reach the surface brightness of ~27 mag arcsec -2 . Thus, we conclude that tidal debris from minor 
mergers can explain the observed dry tidal features in ellipticals at intermediate redshifts, although 
other mechanisms (such as major dry mergers) may also be important. 

Subject headings: galaxies: kinematics and dynamics — galaxies: formation — galaxies: stellar content 
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1. INTRODUCTION 



Recentlv lvan Dokkuml l|2005l hereafter vD05) analyzed 
unprecedentedly deep image of red galaxies at redshift 
around z~ 0.1, and found that they often have diffuse 
tidal features over the scale of ~ 50 kpc. Furthermore, 
there is no associated star formation with these tidal fea- 
tures. vD05 call such features "dry tidal features" . vD05 
also find that 30 % of the red galaxies in their sample are 
undergoing ongoing mergers, and the luminosity differ- 
ence of these merging galaxies are often small (the ratio 
of the luminosities between the bigger and smaller galax- 
ies is more than 0.3). However, it is still an open ques- 
tion whether most dry tidal features are the remnants 
of such major dry mergers or whether other mechanisms 
can produce these features 5 . 

Here, we show that deep images of simulated red galax- 
ies also display red tidal features. The galaxy studied is 
taken from a fully self-consistent A-dominated cold dark 
matter (ACDM) cosmological simulation. We identify 
the responsible incident for the dry tidal feature, and 
find that the feature is due to the tidal tails created from 
an infall of a relatively small satellite galaxy, i.e. a minor 
merger. Our simulation suggests that such minor merg- 
ers may be responsible for some of the dry tidal features 
observed in vD05. 

Our numerical simulations and the simulated ellipti- 
cal galaxy are described in the next section. Section 
presents a quantitative comparison between the dry tidal 
features in the simulated galaxy and those observed in 
vD05, and shows that the simulated features are consis- 
tent with the observed ones. The implications of this 
finding and our conclusions are given in the final section. 
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2. TARGET SIMULATED ELLIPTICAL GALAXY 

The simulated images considered he re are based on the 
ellipti cal galaxy model 2 presented in lKawata fc Gibsonl 
IpM hereafter KG05). KG05 carried out ACDM cos- 
mological simul ations using the galactic c hemodynam- 
ics code, GCD+ l|Kawata fc Gibsonl l2003ajL GCD+ is a 
three-dimensional tree A-body /smoothed particle hy- 
drodynamics (SPH) code which incorporates self-gravity, 
hydrodynamics, radiative cooling, star formation, super- 
novae (SNe) feedback, and metal enrichment. GCD+ takes 
account of the chemical enrichment by both Type II 
(SNe II) and Type la (SNe la) SNe, mass-loss from inter- 
mediate mass stars, and follows the chemical enrichment 
history of both the stellar and gas components of the 
system. 

The cosmological simulation adopts a ACDM cos- 
mology (Q o =0.3, A =0.7, n h =0M9h~ 2 , h=0.7, and 
(T8=0.9), and uses a multi-resolution technique to achieve 
high-resolution in the regions of interest, including the 
tidal forces from neighboring large-scale structures. The 
initial conditions for the simulation s are constructed us - 
ing the public software GRAFIC2 ijBertschingerl 12001^) . 
Gas dynamics and star formation are included only 
within the relevant high-resolution region (~12 Mpc at 
z=0); the surrounding low- resolution region (~43 Mpc) 
contributes to the high-resolution region only through 
gravity. Consequently, the initial condition consists of 
a total of 190093 dark matter particles and 134336 gas 
particles. The mass and softening lengths of individual 
gas (dark matter) particles in the high-resolution region 
are 5.86 x 10 7 (3.95 x 10 s ) M© and 2.27 (4.29) kpc, re- 
spectively. 

The elliptical galaxy found by KG05 (see also 
iKawata fc Gibsonl l2003br) has a total virial mass of 
2 x 10 13 Mq. The galaxy is relatively isolated, with only 
a few low-mass satellites r emaining at z = 0. Figure 1 of 
IKawata fe; Gibsonl (2003b) shows the morphological evo- 
lution of dark matter in the simulation volume and the 
evolution of the stellar component of the target galaxy. 
The galaxy forms through conventional hierarchical clus- 
tering between redshifts z—S and z—l. The morphology 
has not changed dramatically since z=\. 
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KG05 introduced an active galactic nuclei (AGN) 
heating model in the simulation, and found that self- 
regulated activity of the AGN can suppress significant 
late-time star formation - characteristics not encoun- 
tered in traditional dynamical models of ellipticals. As a 
result, their model (model 2 of KG05) succeeds in repro- 
ducing both the observed X-ray and optical properties of 
a massive elliptical galaxy. 

In the chemodynamical simulation, the simulated star 
particles each carry their own age and metallicity "tag" , 
which enables us to generate an optical-to-near in- 
frared spectral energy distribution for the simulated 
galaxy, when combined with our population synthesis 
code (which itself is based upon the p opulation synthe- 
sis models of lKodama fc ArimotdlT99^ . The population 
synthesis models are used to construct the photometric 
images used in our analysis. 

3. DRY TIDAL FEATURES IN THE TARGET GALAXY 

Figures 1 1141 show the morphological evolution of the 
target galaxy over the redshift range z = 0.4 to z = 0.22 
at two different projections: x — y and x — z projection 
in our three dimensional simulation coordinate system. 
The top and middle panels in each figure show the deep 
(R < 27.0 mag arcsec -2 ) and shallow (R < 24.5 mag 
arcsec -2 ) images, respectively. The deep images reveal 
some small tidal features which are not visible in the 
shallow images. The tidal features are similar to some 
of the diffuse tidal tail like features which are observed 
in vD05, e.g., 5-994, 9-360, 10-232, 16-650, 16-1302, 17- 
2819, 21-523, 25-3572 in his sample. For example, in 
Figure ^ the panels at z = 0.37, 0.34, and 0.32 show 
faint features like tidal tails. In fact, these are the tidal 
tails of a small satellite infalling into the target galaxy. 
The infalling satellite galaxy is highlighted by the circle 
in the top-left panel of Figures and [2J The bottom 
panels show the distribution of the star particles. The 
grey dots show the whole particles, while the black dots 
represent the particles originated from the infalling satel- 
lite. The satellite particles arc defined as the particles 
whose radius from the center of the satellite is less than 
20 kpc at z = 0.45, well before the satellite falls into the 
target galaxy. The figures show that the tidal features 
seen in the deep i?-band image correspond to the tidal 
tails of the infalling satellite. The tidal tails of the satel- 
lite start to be developed at z = 0.40. The satellite is 
tidally disrupted at z = 0.37, and the tidal features dis- 
appear around z — 0.25, which corresponds to about 1.3 
Gyr. During the period when the tidal tails are visible 
in the bottom panels, the tidal features also appear in 
the deep image, although, not surprisingly, the strength 
of the features depends on the projection angle. vD05 
claims that the tidal features appear smooth, showing 
little or no evidence for clumps and condensations. On 
the other hand, the simulated tidal features suggest the 
appearance of small clump prior to redshift z=0.34, when 
the satellite has yet to fully disrupt. However, one can 
appreciate from Figure 11 of vD05 that some fraction 
of the observed galaxies present tidal tail-like features 
accompanied with small clumps, for example galaxies 9- 
360, 10-232, 17-2819, and 21-523. They may be an early 
stage of a minor merger. After redshift z=0.34, there are 
appear to be no discernible condensations in the tidal 
tails. 



Another important aspect reported in vD05 is that 
both the host galaxy and the tidal features are red, and 
there is no sign of star formation, i.e., they are "dry". 
Similarly, in our simulation the produced features are 
in fact "dry tidal features" . Figure 03 shows the history 
of the star formation rate of the system, replotted from 
Figure 3 of KG05. Star formation stops in the target 
galaxy just before z ~ 0.4, due to the AGN heating (see 
KG05 for details). The AGN heating also suppresses star 
formation during the infall of the satellite galaxy. 

Tabled displays the luminosity, colors, and mass of the 
target galaxy and the infalling satellite at z = 0.4, i.e., 
before they start merging. Here, we adopt the diameter 
of D = 5 arcsec for a comparison with vD05. The stellar 
mass ratio between the satellite galaxy and the target 
galaxy is about 0.03, which indicates that the infalling 
satellite is small, i.e., a minor merger. This indicates 
that dry tidal features can result from minor mergers. 

Table gives two values for the luminosity of the tar- 
get galaxy: One is the total luminosity, while the sec- 
ond value is the luminosity when the contribution from 
stars which are younger than 3 Gyr is ignored. Since 
star formation stops just before z = 0.4 in the target 
galaxy, the luminosity at z — 0.4 is affected by some 
young stars. As argued in KG05, this late cessation of 
star formation in the simulated galaxy is due to our sim- 
ple modeling of the AGN heating. In this simulation, 
the AGN heating was started at z = 1, which is likely 
too late. If the AGN heating is added at a higher red- 
shift, star formation is expected to stop earlier. For real 
elliptical galaxies, the epoch of the termination of star 
formation is expected to be rather random, and objects 
detected as red elliptical galaxies a re likely at a stage well 
after star formation ceases (e.g., iShiova fc Bekkill998t 
Ivan Dokkum fc Fra nx 2001). Thus, the luminosity when 
the young stellar population is ignored is more appropri- 
ate for objects selected to be red and this is the luminos- 
ity we consider here. For such an object, the luminos- 
ity difference, L x /L 2 = lO^. 1- *^. 2 )/ 2 ' 5 , between the 
target galaxy and the satellite galaxy is 0.05, which is 
consistent with the lowest luminosity ratios found in the 
sample of ongoing mergers in vD05. In fact, those galax- 
ies, 1-2874 and 5-2345 in vD05, show clear tidal tails, like 
those seen at z = 0.37 in Figure ^ Table also shows 
the colors for both the target galaxy and the satellite, 
which indicates that they are both red, and the differ- 
ence in B — R is 0.14. This color difference is driven by 
the difference in metallicity. 

Figure shows the evolution of the luminosity and col- 
ors of the target galaxies. The filled circles indicate the 
colors and magnitudes before z = 0.2, when we found 
the dry tidal features. The open circles represent the 
evolution at later stages. The area enclosed by the dot- 
ted lines indicates the selection criteria in vD05. As dis- 
cussed in vD05, their selection criteria preferably selects 
red galaxies at z < 0.2. The epoch we focus on here is 
a little earlier, and the target galaxy is redder than the 
selection criterion adopted by vD05. However, the epoch 
of minor mergers is random, and such minor mergers are 
expected to happen often in a ACDM cosmology. Hence, 
the results presented here indicate that if minor mergers 
happen at z < 0.2 for elliptical galaxies, the dry tidal 
features observed in vD05 can be produced. 
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Fig. 1. — Deep (upper panels) and shallow (middle panels) predicted _R-band images of the target galaxy at the x — y projections over 
the redshift range z = 0.4 to z = 0.32. Bottom panels represent the distributions of star particles (gray dots), the particles within the 
infalling satellite (see text for details) are highlighted by black dots. The infalling satellite galaxy is highlighted by the circle in the top-left 
panel. The prominent features are indicated by arrows. Numbers shown in the top panels are the tidal parameter, t, (see text) for each 
image. 



TABLE 1 

The luminosity, color, and mass of the target galaxy and the infalling satellite at z = 0.4 a . 





the target galaxy 








the satellite 




Mr 


M fl (age>3 Gyr) b B 


- R c 


M* 


Mr 


B-R 


M* 


(mag) 


(mag) 






(mag) 




(Mq) 


-23.25 


-22.50 


2.53 


3.1 x 10 11 


-19.17 


2.39 


1.0 x 10 10 



a The luminosity is measured within the projected diameter of 5 arcsec, i.e. 26.7 kpc at z = 0.4, while the mass is measured within the 
three dimcntional diameter of 26.7 kpc. 
b The i?,-band luminosity, when the contribution from the stars younger than 3 Gyr is ignored. 
c The color, when the contribution from the stars younger than 3 Gyr is ignored. 



Finally, we quantitatively confirm that the tidal fea- 
tures shown here are consistent with the features ob- 
served in vD05. To this end, we have analyzed the tidal 
parameter, t, which is defined in Section 5.3 in vD05, for 
all the i?-band images presented in Figures 1 1141 Here, 
we briefly explain our procedures to obtain the tidal pa- 
rameter. We refer to the simulated galaxy R-band image 
as "G". First, the galaxy images are fit by an elliptical 
model using the ellipse task in IRAF. The center po- 
sition, ellipticity, and position angle are allowed to vary 
with radius. In cases where there are galaxies which over- 
lap with the target galaxy in projection, the overlapping 
objects are also fit. The final model image is indicated 
by "M". Finally, a fractional distortion image, "F", is 
created by dividing G by M. When the pixel value of F 



is defined by F(x,y), and the mean value of the image 
is F, the tidal parameter t which describes the level of 
distortion is defined as t = F(x, y) — F. Figured demon- 
strates this procedure. 

The values of t are shown in all the top panels in Fig- 
ures We can see that the i?-band images which have 
the tidal features show higher t value (t > 0.1). During 
the epoch when the tidal features of the satellite galaxy 
are visible in the bottom panels (z > 0.25), t appears to 
be higher. As expected, the t values are lower after the 
tidal features disappear at z ~ 0.25. We also find that 
the value of t depends on projection. In some cases, the 
different projections give vastly different t values (com- 
pare for example the x — y and x — z panels at z — 0.32). 
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Fig. 2. — Same as Figure^ but for the x — z projections. The high value of the tidal parameter at z = 0.4 is caused by poor model 
subtraction due to the small separation between the target galaxy and the infalling satellite. 
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Fig. 5. — Time-variation of the star formation rate for the target 
galaxy. 

t = 0.19) disturbed galaxies classified in vD05. Thus, 
quantitatively we confirm that the dry tidal features seen 
in our simulated galaxies are consistent with the dry tidal 
features seen in vD05. 

4. DISCUSSION AND CONCLUSIONS 

We have demonstrated that a minor merger of a 
satellite galaxy whose stellar mass is 10 10 M© with a 
massive elliptical galaxies can produce dry tidal fea- 
tures. In the present simulation, the features last 
about 1.3 Gyr, although the lifetime of such features 
must depend on the orbit of the satellite galaxy. Tra- 
ditionally cosmological numerical simulations of ellip- 
tical galaxy formation have suffered from an inabil- 
ity to su ppress late-time star formation in the simula- 
tion fe.g.. lSuginohara fc Ostrikerll998tlLewis et al.l2000t 




R 

Fig. 6. — Time variation of B — R color and the ii-band mag- 
nitude for the target galaxy. The numbers indicate redshift. The 
filled circles are the colors and magnitudes before z = 0.2, which we 
focus on in this paper. The open circles present the evolution after 
the stages which are discussed in this paper. The area enclosed by 
the dotted lines are the selection criteria in vD05. 



I Kav et al.ll200l iTornatore et alJl2003 IMeza et al.H2003t 
iKawata fc Gibson! I2003bt iBorgani et all l2004j) ~ erebv 
being inconsistent with the observed red color in bright 
elliptical galaxies. However, self-regulated AGN heating, 
of the sort introduced in KG05 aids in the suppression 
of late-time star formation. As a result, both the host 
galaxy and tidal features in our simulation have the red 
colors consistent with those observed by vD05. 

In a ACDM cosmology, minor mergers as seen in 
our simulations are expected to happen much more fre- 
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Fig. 7. — Demonstration of the process to derive the value of the tidal parameter, t (see text for the details). Left panel shows the 
ij-band image; middle panel is the elliptical model; right panel presents the distortion image. 



quently than major mergers. This might explain the 
high detection rate of tidal features in vD05. In addi- 
tion, minor mergers are consistent with the red features 
observed, because mergers with a small galaxy have a 
relatively weak effect on the inter-stellar medium (ISM) 
of the host galaxy. On the other hand, major mergers 
must have a stronger effect on the ISM even for elliptical 
galaxies, where the system is likely to suffer from radia- 
tive cooling. This in turn likely induces star formation, 
resulting in blue galaxy colors. However, some dry tidal 
features should be produced from major mergers, as long 
as these mergers are also "dry". In fact, ongoing dry 
major mergers are observed at various redshifts below 
z=l (e.g., vd05; iTran eFaTIElM iRell et al.ll200fil) . and 
it is known th at such major dry m e rgers can create tidal 
features (e.g.. iCombes et al.lll995l Bell et al.ll2006ft . In 
addition, such major dry mergers do n ot violate the scal- 
ing relations of elliptical ga l axies (e . g . , iNaab et alJl200fit 
iBovlan-Kolchin et al.ll2005l 12006(1 . Our results should 
not be interpreted as degrading the importance of major 
dry mergers, but instead that dry tidal features do not 
necessarily derive from a major dry merger and that a vi- 
able alternate pathway for their appearance may be that 
of a minor dry merger. In other words, infalling galaxies 
over a wide range of mass may be responsible for produc- 
ing the observed dry tidal features in elliptical galaxies. 
Therefore, considerable caution is needed to use the dry 
tidal features for estimates of the major dry merger rate. 
Note that vD05 carefully estimates the major dry merger 
rate, taking into account only the sample of ongoing dry 
major mergers. 

I.Iohnston et al.l l|200lD demonstrate that if the infalling 
galaxy was a dwarf spheroidal with the mass of ~ 10 s 
Mq, like seen in the Local Group, the surface bright- 
ness of the tidal debris would be lower than ~ 28 
mag arcsec~ 2 after the dwarf is well-disrupted. Thus, 
mergers with the most common low mass galaxies are 
apparently not sufficient to produce the features found 
by vD05. Since the infalling galaxy in our simulation 
has a stellar mass of 10 10 M , the surface brightness of 
the resulting tidal debris is much higher, making the fea- 
tures detectable at the depths achieved in vD05. Hence, 
the dry tidal features observed in vD05 are expected to 
be end-products of mergers with a system whose stellar 
m ass is at least ~ 10 10 M ^. Systematic studies like that 
of l.Iohnston et all l)2001j) . but with higher mass satellites, 



would be useful for determining the mass range required 
to reproduce the observed dry tidal features. 

The mechanism which creates the dry tidal feature 
in our simulation is simil ar to wha t is su ggested for 
the formation o f ripp les l|Schweizerl 11980(1 and shells 
l|Malin fc Carterl 119801) observed in nearby elliptical 
galax ies (e.g.. ISchweizer fc Seitzerl Il992t iColbert et alJ 
|2001|). Previous studies demons trate that both a major 
(e.g.. E£incmi&tJL-Sj3er^^[1993) and minor fe.g..lQuin: 
| 1984UDunraz fc CombeslllQSfiUHernouist fc Quinnlll98 
iKoiima fc Noguchil997l) merger can create shell features 
around giant elliptical galaxies. vD05 claims that their 
dry tidal features are different from the ripples and shells 
seen in nearby elliptical galaxies, because the dry tidal 
features look smoother than these other structures. How- 
ever, vD05 also notes that the resolution of his images 
make it difficult to reach strong conclusions about the 
differences between these various types of features. Ac- 
cording t o the previous systematic st udies of minor merg- 
ers (e.g., iHernauist fc Quinnl IT98 8). smoother features 
are produced if the infalling galaxy is more of a veloc- 
ity dispersion supported system, i.e., elliptical galaxies. 
In addition, the tidal features are sensitive to the in- 
fall orbit, and the previous simulations demonstrate that 
off-center infall can create large scale features around the 
host galaxy The infalling galaxy seen in our simulation 
has an off-center orbit, and thus creates the tidal features 
over a large scale, ~50 kpc in radius. 

The dry major merger is considered to be an impor- 
tant mechanism because one major merger is enough to 
explain the observed factor of two increa se in stellar mass 
in bright red galax ies since redshift z=l (Bell et a,Ul2004t 
iFaber et al1 l2005). It is an important question whether 
or not minor mergers, such as those described here, can 
similarly contribute to the increase in the mass of bright 
red galaxies. To check this, we compared the stellar mass 
of the simulated galaxy within a radius of 100 kpc at z=l 
and z=0. As mentioned above, the simulated galaxy suf- 
fers from an excess of ongoing late-time star formation 
(prior to z ~ 0.4 - see also Fig. |SJ). Hence, we exclude 
stars whose age at the present epoch is less than 8 Gyrs. 
As a result, we find little increase in stellar mass subse- 
quent to z = 1 (less than 5% growth). This is because the 
simulated galaxy does not undergo a major merger sub- 
sequent to redshift z = 1.5, and the disrupted satellites 
are much less massive than the host galaxy. (Note that in 
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Table^we measure the stellar mass within a small aper- 
ture and therefore underestimate the total stellar mass 
of the host galaxy.) Therefore, minor mergers such as 
those described here are unlikely to contribute to the ob- 
served mass evolution of red galaxies. This means that 
major dry mergers must occur for a significant fraction of 
the red galaxy population, although we emphasise again 
that while potentially the dominant evolutionary path- 
way, it is not a m utually exclusive one. Interestingly, 
iKong et al.l (|2006) claim that the stellar mass increase is 
imp ortant for the relatively faint early- type galaxies (see 
also lR,enzinill2nn^ . There may be a mass dependence of 
the dry major merger rate. 

Our simulations of a single elliptical galaxy leaves 
many open questions. For example, how often do mi- 
nor mergers happen for bright elliptical galaxies and how 
do the resulting features depend on the mass and mor- 



phology of the infalling galaxy. To address these ques- 
tions, simulations of a large sample of ellipticals will be 
required. 
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